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ABSTRACT 

The investigation of boundary breather states of the sinh-Gordon model restricted 
to a half-hne is revisited. Properties of the classical boundary breathers for the two- 
parameter family of integrable boundary conditions are reviewed and extended. The 
energy spectrum of the quantized boundary states is computed, firstly by using a 
bootstrap technique and, subsequently using a WKB approximation. Requiring that 
the two descriptions of the spectrum agree with one another allows a determination 
of the relationship between the boundary parameters, the bulk coupling constant, 
and the two parameters appearing in the reflection factor describing the scattering 
of the sinh-Gordon particle from the boundary. These calculations had been per- 
formed previously for the case in which the boundary conditions preserve the bulk 
Z2 symmetry of the model. The significantly more difficult case of general boundary 
conditions which violate the bulk symmetry is treated in this article. The results 
clarify the weak-strong coupling duality of the sinh-Gordon model with integrable 
boundary conditions. 
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1 Introduction 



In the study of the sinh-Gordon model restricted to a half-hne by integrable boundary conditions, 
an important issue is to determine how the reflection factor depends upon the two parameters 
introduced at the boundary. Some progress was recently made in tackling this question ||, 
^. Namely, the idea of was to compute the bound-state spectrum of the model in two 
different ways, firstly by using a boundary bootstrap principle, and secondly by quantizing the 
classical boundary breather states using a WKB approach |]^, [T^. Each method provides an 
independent description of the energy spectrum of the boundary states, and comparing the two 
provides information relating the boundary parameters and the bulk coupling to the reflection 
factor. Clearly, one expects the most general situation where the two boundary parameters are 
independent to be technically more involved. It is therefore natural - as was done in - to first 
implement the programme outlined above in the special case where the bulk symmetry — > — 
is preserved at the boundary, by requiring the two parameters to be equal. A further step was 
taken in where a perturbative calculation to lowest order in the bulk coupling and to first 
order in the difference of the two boundary parameters was used to make an informed guess as 
to the general dependence of the reflection factor on those boundary parameters. The results 
of the present paper, which extends the analysis in to the case where the two boundary 
parameters are different, will underpin that guess. Finally, a version of a weak-strong coupling 
duality transformation appropriate to the model with boundary conditions was also proposed 
in [0] and the conclusions drawn from our paper give further support to that conjecture. 

Following a brief introduction to the sinh-Gordon model with integrable boundary conditions 
in section 1, we recall the classical breather solutions in section 2 and describe the properties 
and facts we need for the semi-classical quantization. In section 3, we review the boundary state 
bootstrap and use it to derive a formula for the boundary state energy spectrum. The adapted 
Dashen-Hasslacher-Neveu method is described in section 4 and used to provide an alternative 
calculation of the boundary state spectrum in terms of the boundary parameters. Our results 
and some additional remarks are summarized in section 5. 



2 The sinh-Gordon model on the half-hne 

The sinh-Gordon model describes a single real scalar field in 1+1 dimensions with exponential 
self-interaction. The field equation is 

- dl4> + sinh(v^/30) = 0, (2.1) 

where m and /3 are parameters and we have used normalizations customary in affine Toda field 
theories of which the sinh-Gordon model is the simplest 0. The dimensional mass parameter 
m will be set to unity for convenience. 
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In contrast to the sine-Gordon model, with its sohton and breather solutions, the sinh- 
Gordon model is at first sight relatively uninteresting. There is a constant vacuum solution 
= and, in the quantum theory, the small oscillations around this vacuum correspond to 
the sinh-Gordon particle. In the bulk, the spectrum of the model consists of a single species of 
scalar particle interacting with itself. Nevertheless, however simple the model may appear, it is 
not easy to analyze it directly |T^, and much of what is known about it has been deduced from 



features of the lightest breather in the sine-Gordon theory. 

The sinh-Gordon model is integrable which implies in particular that there are infinitely 
many mutually commuting, independent conserved charges Q±s, where s is any odd integer, 
and the S-matrix describing the scattering of two sinh-Gordon particles with relative rapidity 



9 is conjectured to be given by [T^, |T9 



In (U-'4) we have used the convenient block notation 



sinh ( - + 

(3;) = , (2.3) 

^ ^® sinh(f-^)' ^ ^ 

and the coupling constant B is related to the bare coupling constant (3 hj B = 2/5^/ (47r + (3^). 

For compactness, we will generally omit the subscript O from the block notation since it is 
generally clear from the context what is intended. 

The sinh-Gordon model can be restricted to the left half-line —00 < x < without losing 
integrability by imposing the boundary condition 

5x0lo = — ^ (^eoe ^2"^' ' ' - eie^"^' ' ' J , (2.4) 



where Eq and ei are two additional parameters |12, Again, m sets the scale but is taken to 



be unity for convenience in what follows. This set of boundary conditions generally breaks the 
reflection symmetry — —0 of the model although the symmetry is preserved when Eq = Si = e. 

Assuming factorization, the description of the sinh-Gordon particles on the half line does not 
only require the two-particle scattering amplitude ( ^.2| ), but also the amplitude for the reflection 
of a single particle from the boundary. This reflection amplitude was deduced from the lowest 
breather reflection amplitude in the sine-Gordon model by analytic continuation in the coupling 
constant (i.e., the continuation A —2/B in the notation of |]12|)- Using the breather reflection 



amplitudes calculated by Ghoshal [|13l, the analytic continuation leads to | 

K(ee e 3) = (1)(2 - i?/2)(l + i?/2) 

' °' ''^^ (l-E(£o,^i, /?))(! + i?(^o,^i,/3))(l-F(£o,^i,/3))(l + F(eo,ei,/5))' ^ ' ^ 



where we are again using the block notation from (|2.3|) but, here, 9 represents the rapidity of a 
single particle. When the bulk reflection symmetry is preserved one of the two parameters E or 

'^In Ghoshal's notation E = Brj/n, F = iBd/iT. 
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F vanishes. Without loss of generahty, we shall take the vanishing parameter in the symmetric 
case to be F. All reflection factors satisfy the crossing-unitarity relation which, in the case of 
scalar reflection factors, reads, 

[e + K, [e - '-^^ 3(26) = 1. (2.6) 

Actually, ( p.5| ) is the simplest solution to the crossing-unitarity relation using the S-matrix (p.2|), 
taking into account the independently calculated classical limit of the reflection factors. In |^ 
the classical reflection factor was found to be given by the formula 

K^ie, So, £i, /3) = - ^ — ^TTT V (2-7) 

(1 — ao — ai){l + flo + — flo + + ao — ai) 

in which it was convenient to use an alternative expression for the boundary parameters, viz., 

eo = cosvrao, ei = cosvrai. (2.8) 

The formula is well-defined provided we select the ranges < < 1 for i = 0, 1, to ensure a one 
to one correspondence between the alternative parameters. Clearly, ( p.5| ) has the correct limit 
provided 

E—^ao + ai, F ao — fli. (2.9) 

Recently, on the basis of a perturbative calculation, the relationship between the various pa- 
rameters was conjectured to be [^], 

E = (ao + ai)(l -fi/2), F= (ao-ai)(l-5/2). (2.10) 

The principal purpose of this paper is to provide further evidence for these formulae. 

It was noted in that contrary to the situation on the whole line, the sinh-Gordon equa- 
tion restricted to a half-line by integrable boundary conditions has non-singular, finite energy, 
breather solutions. These solutions were described in some detail in 0, particularly for the 
special case in which eo = ci. Here, we shall concentrate on the solutions containing two in- 
dependent boundary parameters. The existence of these special solutions and their associated 
states makes it clear that the sinh-Gordon model has a rather more interesting structure than 
one would be led to believe on the basis of bulk calculations. 



3 Boundary Breathers 



The boundary breathers may be conveniently described following Hirota's prescription ||l4 
Generally, provided we set 

(3.1) 
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and choose the two r-functions as follows: 



r± = 1 ± (El + ^2 + E3) + (^12^1^2 + ^13^1^3 + ^13^1^3) ± A^2ArsA23E^E2Es, (3.2) 
with 

Ep = e'^f^+^*+'=p, flp = 2 coshpp, bp = 2 sinhpp, Apg = tanh^ (^^^y^) ' ^^-^^ 



then we have just enough freedom to accommodate the general boundary conditions ( p.4|) . In 
effect, with general boundary conditions, the solutions have the flavour of a 'breather' superim- 
posed on a stationary 'soliton', borrowing the language of the sine-Gordon model. In detail, it 
is enough to choose Ei and E2 to be complex conjugate partners and periodic in t, and to take 
E^ to be real and independent of t; thus, setting pi = — p2 = ip and ps = 0, one has: 

^2x cos p+2it sin p+c ^ ^.^ = ^2x+d. = - tau^ p, A13 = A23 = - tau^ (p/2) . (3.4) 

The period of the breather is tt/ sin p. To match the boundary conditions, c and d need to be 
determined. It is convenient (actually just a change of origin in t) to take c to be real. Then, 
the boundary conditions ( p^ ) are satisfied provided 



tanp' tan^(p/2) 

where 



e'^ = (3.5) 



sin - sin ^ 
sin ^ + sin ^ 



2 1 + COS p COS ^ + cos p cos ^ " — COS p 

~l-cosp COS - cosp COS "("°-"^) + cosp" 



(3.6) 



In ( p. 61) we have made use of the alternative parametrisation (p.8|) . There are other ways of 
writing ( |3.6| ) which are useful when it comes to evaluating certain integrals. For example, 
setting q = tan(p/2) and q± = tan(7ra±/2), where a± = (qq ± '^i)/2, we have 



= J_ IlL% izi^, . = (37) 



For symmetrical boundary conditions with £0 = ^1 = ^1 the terms containing E^ are not 
required (effectively d —>■ —00), and the other pieces of the r-functions collapse to: 

= 1 ± 2cos(2tsinp)e---^,/i±^ - ( '-±^] . (3.8) 

tan p \ e — cos p V ^ ~ cos p / 
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The expression ( |3.8| ) is quite easy to analyze. In particular, the solutions ought to be real and 
have no singularities in the region x < 0. It was pointed out in that these requirements are 
met provided the parameters e and p are suitably restricted: 

— 1 < £ < and cosp<— e. (3.9) 

It is interesting that the upper limit cos p = —e corresponds to a choice of frequency at which the 
amplitude of the solution has collapsed to zero yielding the vacuum configuration, = 0. This 
is quite unlike the bulk breathers of the sine-Gordon theory which approach the vacuum as their 
frequencies approach zero. On the other hand, this behaviour of the boundary breathers is closer 
to that of a standard harmonic oscillator of a given frequency whose classical amplitude may 
be arbitrarily small. Once quantized, the energy spectrum of the oscillator depends only on its 
frequency and the zero-point energy is there to remind us that an oscillator of arbitrarily small 
amplitude is distinct from the vacuum. Similarly, we expect that the sinh-Gordon boundary 
breathers will have a non-zero ground state energy. 

For general boundary conditions the properties of the solutions are less amenable to analysis. 
However, numerical investigation of the general boundary breathers indicates that they are non- 
singular in the region a; < provided 

7r(ao + ai) 7r(ao + ai) 7r(ao - ai) , . 

cos < 0, 0<cosp<— cos , cos > 0. (3.10) 

Again, the breathers have frequencies bounded below because the parameters are restricted. We 
shall take the parameters to satisfy Oq > ai (since one must be larger than the other unless they 
are equal, we take the larger to be ao), and to lie within the region defined by: 

0<ao,ai<l, 0<a_<^, ^ < a+ < 1, 7r(l - a+) < p < ^. (3.11) 

As we have noted, the boundary breathers for boundary conditions preserving the symmetry of 
the sinh-Gordon equation are included as the special case oq = ai. 

The energy functional of the sinh-Gordon model incorporating the boundary condition (|2.4|) 
is given by 

S[<P] = £ dx Q^^ + ^0'^ + ^ (cosh(v^/?0) - l)^ 

+ ^ (.o(e-^^(°'*^ - 1) + .,(e^*(°'*) - 1)) , (3.12) 
but it is most easily calculated in terms of the r functions as a boundary term , 



Using this, the energies of the boundary breathers were calculated in [B[ and are 



(3.13) 



(-2 - 2 cosp + (sin ^ + sin ^)') . (3.14) 
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For the symmetric case with ao = oi the energy simphfies to 



^breather = -^l- cos p-e). (3.15) 

In the quantum field theory, the continuum of boundary breather solutions is expected to 
lead to a discrete spectrum of boundary bound states. To estimate this spectrum we shall follow 
the method described in 0] which adapts the techniques developed by Dashen, Hasslacher and 
Neveu (DHN) 0. Although it is not obvious that this method gives an exact result, it is 
clear that the result is non-perturbative, in the sense of being an all orders computation in 
perturbation theory in terms of the bulk coupling constant and the boundary parameters. 

One of the ingredients to the DHN prescription is the classical action computed over a single 
period of the boundary breather. This quantity is relatively straightforward to calculate for the 
symmetric boundary condition but, as far as we can see, is not tractable analytically for the 
general breather. Nevertheless, we have conjectured the result and checked it numerically in two 
different ways. The principal difficulty lies in calculating the integral representing the kinetic 
energy over a single period. We maintain this should be given by a simple expression linear in 
a+, that is 

T rO _ 

dt I (ix 0^ = — (p — 7r(l — a+)), (3.16) 







where T = it / sin p is the period of the general breather. When oq = cti, eq (|3.16| ) agrees with 
the result given in @; using Maple it is possible to check that the integral is independent of 
a_ by differentiating the integrand with respect to a_ and calculating the resulting integral 
numerically, obtaining 0; it is possible to check directly, again by numerical integration within 
the ranges of parameters ( p.lO| )), that / depends linearly on both p and a+. Some details of 
these computations will be given in appendix A. The conjectured expression (|3.16|) seems to us 



to be an astonishing result for which we are quite unable to find an analytical derivation, despite 
its beguiling simplicity. 



4 The boundary bootstrap 

For certain ranges of the parameters E and F the particle refiection amplitude (|2.5| ) has simple 
poles at particular imaginary values of 6 on the physical strip, < Im(6') < 7r/2. These are 
expected to be due to the propagation of virtual excited boundary states although there are 
also other potential explanations via generalized Coleman-Thun mechanisms |P, |T^. The reason 



for this is the following. Once the boundary condition breaks the bulk symmetry of the model, 
the lowest energy field configuration is no longer = 0. This means that in the bulk there will 
be effective odd-point couplings in addition to the standard even-point couplings ||^ and these 
may be used to construct loop Feynman diagrams 'attached' to the boundary, some of which 
may generate poles in the reflection factor Kg. This possibility is difficult to analyze and we are 
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unable to pursue it here. We shall simply assume the poles are due to virtual excited boundary- 
bound states without concerning ourselves with their dynamical origin. The amplitudes for the 
reflection of the sinh-Gordon particle from these excited boundary states are obtained by the 
boundary bootstrap |Tl|, ^. When the reflection factor ( p.5|) has a pole ai 6 = iip with 



<■?/'< 7r/2 then the reflection factor corresponding to the associated excited boundary state 
is calculated via the relation 

Ki{d) = KQ{d)S{d-iij)S{d + iij), (4.1) 

where S{9) is the two-particle S-matrix ( p^.2| ) and Kq{9) is the ground state reflection factor. 
Also, since energy is conserved, the energy of the excited boundary state relative to the ground 
state is given by 

^1 = £:o + m(/5)cos^, (4.2) 

where m{(3) is the mass of the sinh-Gordon particle. 

As a parenthetical remark it is worth mentioning the consequences of this bootstrap proce- 
dure for a free, real scalar field of mass m with a linear boundary condition at a; = 0: 

dx<P\Q = -mXcj). 
Its S matrix is unity but the reflection factor is given by 

i sinh 9 + \ 1 



K 



2sinh^-A (l + 2a)(l-2a) 



where A = cos avr in the second formula. This has a pole which indicates (provided — m < A < 
or 1 > a > 1/2) a boundary breather of a fixed frequency tu, given by uj"^ = m? — A^ = 
w? sin^ avr. Explicitly, the appropriate normalizeable solution to the Klein-Gordon equation and 
the boundary condition is = Ae~^^ cosujt. Repeated application of the bootstrap equations 
leads to a tower of boundary states, each with the same reflection factor, but with energies given 
by 

Sn= £-Q + nuj. 

This is just as one would expect for a harmonic oscillator attached to the boundary although an 
alternative dynamical argument would be needed to determine the ground state energy £q = uj /2. 

Returning to the sinh-Gordon case, let us begin by considering the classical reflection factor 
( |2.7|) in the light of the boundary breather parameter restrictions (|3.1CI| ). It is immediately clear 
that because of the restrictions on the parameters only one of the four factors in the denominator 
of ( p.7| ), namely (1 — ao — ai) has a pole in the physical strip, (remember, we have taken the 
principal values 1 > 0+ > 1/2, 1/2 > a_ > 0). We may also recall the perturbative result 
reported in |Q which implies that at least for sufficiently small coupling E and F do not roam 
far from their classical values. Bearing in mind the classical limits ( p.9| ), this remark suggests 
that we may consider the poles associated with E alone, ignoring F. 
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The remainder of the discussion in this section follows closely the calculations reported in 
but, for the sake of completeness and small changes of notation, they will be repeated here. 
The regions in E where the amplitude (|2.5|) has poles on the physical strip are 



I : 2 > E > 1 and II : -2 < ^ < -1; (4.3) 

since < i? < 2, the other factors never have poles in the physical strip. In region I, ip = 
7r{E — l)/2 and, using ( [4.1| ), we derive the reflection factor for the first excited state. 



K,= {1 + E + B)il-E-B) 

' {1-E){1 + E) {1-E + B){1 + E-B)' 

where the 'spectator' factors have been lumped together in 

(l)(l + i?/2)(2-i?/2) 

The reflection factor ( [4.4| ) in turn has a new pole at ip = tt{E — 1 — B)/2, provided B < E — 1, 
indicating another excited state whose reflection factor is 

^ ^ Ks (i + E + B){l-E-B) 

' {1-E + B){1 + E-B) {1-E + 2B){1 + E-2B)' ^ ^ 

Continuing the procedure leads to a set of excited states with associated reflection factors given 

by, 

^ Ks il + E + B)il-E-B) 

" {l-E + {n- 1)B){1 + E-{n- 1)B) {1 - E + nB){l + E - nB)' ^ ' ^ 

Note, the pole corresponding to the (n + l)st state will be within the correct range provided E 
satisfies 2 > ii^ > 1 + nB. Thus, for a given E and B there can be at most a finite number of 
bound states, and possibly none. Note too that the reflection factor for scattering from the nth 
bound state also contains a pole corresponding to the {n — l)st bound state. Since the factor 
Ks is merely a spectator, at no stage will poles with F-dependent positions enter the game if 
they did not do so at the start. 

The energies Sn of the boundary states are found by repeatedly applying ([4.2|). They are 
given by, 

Sn+i = Sn + m{(3)cos^{nB-E + l). (4.8) 

This is the result that we will compare with the quantization of the classical breather spectrum 
in order to determine how E, F and m{(3) depend upon and ai. 

The poles in region II do not represent a new set of states. In the discussion section we shall 
make some comments concerning the relative roles of E and F. 



8 



5 Semi-classical quantization 



The first step in carrying out the semi-classical calculation is to solve the sinh-Gordon equation 
linearised in the presence of the boundary breathers. Setting (j) = (po + rj, where (po is a classical 
breather, the linear wave equations which ought to be satisfied by the fluctuations are: 



^ - + 4r/ cosh = 0, (^^ + v [eoe"^^"/^ + e,e^^^'/^'] ^ = 0. (5.1) 



x=0 

It is convenient to solve (|5.1| ) by perturbing ( |3.1| ). In other words, we may take 

V = , 

r+r_ 



(5.2) 



with St± chosen in turn by adding a pair of 'small' exponential terms to ( p.2| ) as follows. Take 
Ci and 62 defined by 

ei = Ai 6-^"^*+^^^, 62 = A2 e-^'^*-'^^ uj^-k'^ = 4, 

where Ai, A2 are infinitesimally small. We use Hirota's method again, but with five basic ingre- 
dients this time, instead of the previous three, keeping only the terms linear in 61 and 62. Thus, 
we write, 

St± = ^ 6p f ± 1 + {jJplEi + /ip2-E2 + /ipS-E's) 
p=l,2 ^ 

±(/^pl/^p2^12-£^l-E'2 + /^pl/^pS^lS-^l-^S + /^p2/^p3^23-£^2-£^3) 

+/ipl/ip2/ip3^12^13^23-El-E2-E3^ . (5.3) 

The new coefficients obtained from the general formulae ( |3.3| ) and given by: 

1 ik cos p — uj sin p — 2 1 ik cos p + u sin p — 2 

/^ii = — = -, ■■ rr;^ /^i2 = — = ] : (5.4) 

yU22 "ik COS p — ujsmp + 2 p2i ik cos p + ujsm.p + 2 

1 _ik-2 
/i23 tk + 2 

Matching the boundary condition at x = fixes the ratio A2/A1 to be Q 

^ A2 ^ik — 2 cos p\^ ik — 2 ik — 2c+ ik + 2c^ {5 5) 

Ai \ik + 2 cos p J ik + 2 ik + 2c+ i/c — 2c_ ' 



Remark: the right hand side of the corresponding result in eq(5.5), has a misprint and the printed 
formula ought to be inverted to obtain the correct result for Kb 
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where c± = cos7r(ao ± a.i)/2. The result ( |5r5| ) is not intended to be obvious; its derivation rehed 
heavily on using a symbolic algebraic manipulation language — in our case, Maple. In the limit 
X — > —oo, the fluctuation r/ is a superposition of left and right moving plane waves, 

r/ ~ Ai e-'"* (e''^" + Kb e-''^") , (5.6) 

and the relative phase of these waves defines the classical reflection factor corresponding to 
the boundary breather. Taking cosp = — c+, the breather collapses to the static ground state 
solution and the reflection factor collapses to 

^ ik — 2 ik + 2c+ ik + 2c_ 
ik + 2 ik — 2c+ ik — 2c_ ' 

which, with k = 2sinh6', is simply an alternative form (ie without using the 'block' notation) of 
the expression ( p.7|) . 

The period T = vr/ sinp of the boundary breather defines the 'stability angles' via 

ri{t + T,x) = e-'^riit, x) = e-"^^r]{t, x) (5.8) 

and the field theoretical version of the WKB approximation makes use of the stability angles 
together with a regulator to calculate a quantum action. If there are no boundaries, it is natural 
to add some artificially to render the spectrum discrete and facilitate the necessary calculations. 
One way, the simplest and most commonly used, would be to place the field theory in an interval 
[— L, L] with periodic boundary conditions and to manipulate the sum over the discrete stability 
angles. However, since we have one boundary already prescribed, we need to do something else 
as suggested in 0. It is convenient to treat the sinh-Gordon model in the interval [— L, 0] and 
to impose the Dirichlet condition //(t, —L) = 0. Since the limit L —>■ oo will be taken eventually, 
the stability angles for the boundary breather (i^b), or the vacuum solution (z/q) are effectively 



determined by the reflection factors given in ( |5.6|) or (|5.71 ), respectively. A potentially more 



interesting calculation but one which we do not yet have the machinery to carry out, would 
be to consider two sets of two-parameter boundary conditions. However, the breathers we are 
considering are inadequate for that. 

Following H, [l^ we need to calculate a sum over the stability angles and use it to correct 
the classical action. Thus, 

A=^-J2i'^B-'yo)^^J2 (\Al + 4-^^), (5.9) 
where /c^ and ko are the sets of (discrete) solutions to 



e 



2ikBL _ f iks + 2 COS p\ ikB + 2 ikB + 2c+ ikB -2c. 



ikB — 2 cos p J ikB — 2 ikB — 2c+ ikB + 2c. 



^2ikoL + 2 iko- 2c+ iko - 2c. 



^Mk,L ^ _ (5.10) 
ik^ — 2 zkq + 2c+ zfco + 2c_ 
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Once A is known, the quantum action is defined by 

S'q„ = ^ci-A, (5.11) 

where the classical action is readily calculated from the kinetic energy integral ( |3.16| ) and the 
total energy ( p.l4|) : 

Sd= I dt I = ^ ( p - 7r(l - a+) + ( cosp + 1 - i(l - c+)(l + c_) ) ) . 

Jo J-oo P V smp V 2 J J 

(5.12) 

We shall proceed along the lines described in ||^ noting that for large k the solutions to either 
of ( |5.1CI| ) are close to 

1\ TT 

and so it is reasonable to set (kB)n = (^o)n + i^{{kQ)n)/L where, for L large, the function n is 
given approximately by 

^2Mk) ^ f ik + 2 COS p ik + 2c+ Y ^3^, 
\ik — 2 cos p ik — 2c+ / 

Interestingly, all dependence upon c_ has dropped out and therefore, from this point on, the 
calculation is identical to the corresponding part of the calculation presented in |Q with e = 
cosvra replaced by c+. 

In terms of k, the expression (|5.9| ) is rewritten 

A ~ — y (^° )"^((^")" ) + 0(1/l2), 

and this, in turn, as L oo can be converted to a convenient (but actually divergent) integral, 

T /•°° „ kK(k) 
A = — / dk , ' 5.14 

with which we shall have to deal. Note that k vanishes when cosp = —c+. Integrating (|5.14|) 
by parts we find 



T 



2ti 

where 



A = — + 4 -/ dk—VWn\, (5.15) 



dk 



dn _ 4 cosp 4c+ 

dk A;2 + 4cos2p A;2 + 44' ^ ' ' 
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and we note that with a suitable choice of branch 

4 cos p 4c+ 



as k ^ oo. (5-17) 

Using ( p. 171 ) and recalhng that cosp < — c+, we deduce that k approaches zero from above as 
k oo. Also, from ( ^.16|) it is clear that the derivative of k is positive near k = but negative 



as A; ^ oo. Hence, the first term in ([5. 151) is well-defined and the appropriate branch of k, has 



k(0) = 0. On the other hand, the derivative of k is not decaying sufficiently rapidly to ensure the 
second term in ( ^.15| ) is finite. However, this was to be expected since a perturbative analysis 
of the sinh-Gordon model confined to a half-line needs mass and boundary counter terms to 
remove logarithmic divergences (which would be removed automatically by normal-ordering the 
products of fields in the bulk theory). With this in mind, the integral remaining in ( [5. 151 ) should 
be replaced by 

4 cos p 4 cos p 4:C+ 4c-( 



V A;^ + 4 cos^ p A;2 + 4 A;^ + 4c+ A;^ + 4 / 

the first counter-term removing the bulk divergence and the second being there to remove a 
similar divergence associated with the boundary. In effect, we are regarding the parameter a 
as describing the bare coupling which appears in the boundary part of the Lagrangian once it 
is written in terms of normal-ordered products of fields. The counter-terms clearly respect the 
symmetry and the whole expression vanishes when p = 7r(l — a_|_). The integrals in ( [5. 18] ) need 
to be treated carefully noting that cos p > but c+ < 0. 



Besides the towers of real solutions to (|5.10| ), there is also a discrete set of solutions for 



which fco and ks are pure imaginary. These were not discussed in |0 but including them in the 
argument would not have altered the conclusions, as we shall show. First, notice that as L — > cxo 
the imaginary solutions which survive are either zeroes or poles of the right hand sides of the 
equations in (|5.1CI| ), according to the signs of ik^ or iks- Thus for k^ we have ik^ = 2, — 2c4.if 



iko > and iko = — 2,2c+ if iko < 0, while for ks we have ikB = 2 cosp, 2 if ikB > and 
iks = —2 cosp, —2 if ikB < 0. However, the signs should be disregarded because for either sign 
each solution for a specific value of |iA:o| or \ikB\ represents a single 'bound-state' function rj. 
Taking this into account, these special solutions contribute to A an additional piece: 

rrt ■ ■ \ 7rsin7ra+ 

I (sm p — sm 7ra_|_) = vr . (5.19) 

sin p 

Assembling the various components leads to 

2 

A = TT (cosp + cos7ra+ + psinp + 7r(a+ — 1/2) sin7ra+) . (5.20) 

sm p 



Recalling (|5.12| ), and using ( ^.201 ), the quantum action defined in (|5.11|) is given by an expression 



of the form 



4 /cosp 7r\ Stt / 7r\ T(a+,a^) 

S^u = -b(—^ + P-7^)+-^ vra+ - - + ^^t^, 5.21 
±5 \ sm p 2 / \ 2 J sm p 



12 



where T is independent of p, 

r(a+,a_) = 2n ( - c+)(l + c_) + + — + (a+ - 1/2) sinvra^ 
\ p p'^ TC 

Once the quantum action is determined, the quantum energy is defined by 

^ _ agqu _ sin^ p dSqu _ 4 r(a+,a_) 
ai vr cos p op ttB tt 



and the WKB quantization condition states that 

7r\ Svr / vr 
5 2)^]3 



Wc,u = S^, + T£^^ = -(p--) + — (na+--)= 2nn. (5.23) 



Where, n is a positive integer or zero.| Hence, the energies of the quantized boundary breather 
states are determined by a set of special angles pn, 



Pn = |(l+ni?-^(2a+-l)), (5.24) 



and given by 



4 r 4 TT / 2nB ^ x\ r , , 

tn = ^cospn = FT COS - [ uB + 1 - 1) • (5.25) 

ttB tt tvB 2 \ p^ / ^ 

Notice that as /? — 0, p„ ^ 7r(l — a+) independently of n. Thus, the frequencies collapse to 
the lowest allowed frequency, namely cuq = 2sin7ra+. On the other hand, in the same limit the 
energies are independent of /3 and non-zero, 

Sn-^[n + oJo. (5.26) 
This is precisely the spectrum of a harmonic oscillator vibrating at the fundamental frequency 

UJQ. 

Using ( [5.25| ) the corresponding differences in the energy levels are given by 

£n+i = ^„ + — sin — COS - (^^(2a+ - 1) - (^n + - j Sj . (5.27) 

Comparing ( |5.27| ) with the outcome of the bootstrap calculation ( [4.8| ) allows us to identify the 
parameter E which appeared in the expression for the refiection factor (|2.5| ). Thus, from the 
first excited level we deduce, 

E(ao, ai, (3) = ^(2a+ - 1) - ^5 + 1 = 2a+ (^1 - |^ = (ao + ai) (l - f ) . (5-28) 



Including the imaginary solutions of (5.10) obviates the need for making the change n — > n + 1/2, since we 
shall see below the zero-point energy is automatically correct. In this sense, including the imaginary solutions 
leads more naturally to the expected result. 



13 



and then all the other levels match up without further restriction. This is in agreement with 
the suggestions made in 0. 

As was pointed out previously in the comparison with ( [4 .81 ) also permits us to deduce 
an expression for m{j3), the mass of the sinh-Gordon particle: 

m(/3) = — sin— . (5.29) 



6 Discussion 



From the point of view of the classical reflection factor ( |2.7|) , it is clear that the parametrisation 
( |2.8| ) is the most appropriate. Moreover, the expressions for Eq and Ei in terms of Oq and Oi are 
invariant under independent changes of signs of either oq or ai and this symmetry is incorporated 
in the factors appearing in the denominator of ( p.7|) . We expect the symmetry under reversing 
the signs of oq and ai should persist in the quantum reflection factor, and indeed that was part 
of the thinking behind the conjectured forms for E and F presented in ( p.lO|) . In other words, 
having calculated E in ( |5.28| ) we can immediately deduce the partner expression for F. 

We found that the general breathers exist as non-singular real solutions provided the param- 
eters satisfy the constraints ( |3.11| ). It is gratifying to discover at the end of the calculation that 
the renormalization factor 1 — B/2 in (|2.10|) for real coupling f3 is restricted to lie in the range 

B 

< 1 < 1, 

2 - ' 

and thus has the effect of scaling the parameters but preserving the constraints ( |3.11|) . Notice 
that with E and F given by (|2.1CI|) it is impossible to have E and F lying simultaneously in 
regions ( [4.3|) for which the bound state poles of the reflection factor lie on the physical strip. 
Thus, the poles depending on or F are alternatives and actually lead to the same tower of 
states simply seen differently in terms of the parameters ao and ai. 

Changing the sign of either oq or ai takes us outside the principal ranges of these parameters 
and obliges us to reformulate the breather solutions. In fact, there are no others besides the 
ones we have found although the expressions for them will be a little different if we adopt a 
different principal region. 

It is worth emphasizing that the expressions (|2.1(]|) incorporate a weak-strong coupling du- 
ality 0] which extends that enjoyed by the S-matrix itself. If a new triple of coupling constants 
{al,al,l3*) is defined by 

Air 

{al,al,p*) = — (ao,ai,/5). 



then 

E(nl n* R*\= ^^^^^ , ^^^^ 



47r B 

E{al, al, (3*) = -^{ao + a^)- = (ao + ai)(l - B/2) = E{ao, ai, /?), (6.1) 
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and similarly for F. 

We must acknowledge that the semi-classical approximation, although non-perturbative, is 
not guaranteed to be exact. However, for the bulk sine-Gordon model, the Dashen, Hasslacher 
and Neveu approach does give exact information concerning the spectrum and we would ex- 
pect that the same should be true here. Unfortunately, alternative exact computations of the 
spectrum are not yet available for comparisons to be made. 

As far as other models are concerned, it will be interesting to try out these methods in other 
cases where the boundary parameters form a continuous set. The simplest such example, in 
which there are two distinct one-parameter families of boundary parameters, is the model based 
on the a^^ root data 
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A Appendix 

In this appendix we shall give some information concerning the kinetic energy integral (|3.16|) . 
The first step in evaluating the integral is to perform the integration over time. This is relatively 
easy and leads to: 

^ 2 _ Svrsinp \2CVB'^ - A'^ ^ B _ 2A^JD'^ ~ ^ 



I{x) = / dt^ 







AD-BC ^B'^-A^ AD-BC ^D'^-C'^ 



(A.l) 

where the right hand side is built from the components A, i?, C, D which are defined by 

A = —(1 - g2)(i _ rv), C = - — (1 - g2)(l + rv), 
q q 

B = — ((rv — s^q^u^) + q^(l — s^q^rvu^)) , 

D = ^ (~(rv + s^q\^) + q^(l + s^q\vu^)) , (A.2) 

q2 \ I 
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with 



u = e'^=°'^^ V = e'^. (A.3) 

The parameters q,r,s were defined previously in ( p.7| ). Notice, when = cii we have r = 
allowing all of these expressions to simplify dramatically. Consequently, the integral of the right 

hand side of ( |A.1| ) can be done by making a suitable change of variables. Thus, for 
we find: 



dxl{x) = —{p-n{l- a)). (A.4) 



2 

oo 



However, in the general case we have not found a change of variables which simplifies ( [A.l| ). 
Therefore, we have had to proceed numerically. 

One check we have made is to differentiate ( |A.1| ) with respect to a_ and demonstrate 



_d_ 

(9a_ 





dxl{x)=0, (A.5) 

-oo 



provided a± lie within the ranges ( |3.11 ). Another is to calculate the integral directly as a function 



of a+ and p and demonstrate it depends linearly on each of these parameters separately. The 
following plots generated by Maple indicate unambiguously the linear dependence on a+ (Figure 

1) 




Figure 1: The integral as a function of a+. 

or on p (Figure 2) 
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Figure 2: The integral as a function of p. 

In Fig(l) we have overlaid several plots. The dots represent the numerical values of the integral 
as a+ varies from 1 — p/vr to 1 — a_, in accordance with ( |3.11|) , for each of five different values 
of p (viz. p/t^ = (2 + A;)/16, k = 1, . . . , 5); the continuous line is simply the conjectured value 
of the integral given in ( |3.16| ) plotted as a function of a_|_ for the same five values of p. In all 
plots, the factor Svr//^^ has been ignored. The vertical line indicates the upper bound on a+ for 
a_ = 1/8. Fig(2) presents similar plots for the integral as a function of p/vr in the range 1 — a+ 
to 1/2 for five different values of a+ (viz. a+ = (8 + /c)/16, /c = 1, . . . , 5). We regard the sets of 
parallel lines depicted in Figs(l,2) as very convincing numerical evidence for ( [3.161 ). 
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